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We review the theory of second-order (ferro-) elastic phase transitions, where 
the order parameter consists of a certain linear combination of strain tensor 
components, and the accompanying soft mode is an acoustic phonon. In three- 
dimensional crystals, the softening can occur in one- or two-dimensional soft 
sectors. The ensuing anisotropy reduces the effect of fluctuations, rendering the 
critical behaviour of these systems classical for a one-dimensional soft sector, and 
classical with logarithmic corrections in case of a two-dimensional soft sector. The 
dynamical critical exponent is z = 2, and as a consequence the sound velocity 
vanishes as c s oc \T — T c \ 1 / 2 , while the phonon damping coefficient is essentially 
temperature-independent. Even if the elastic phase transition is driven by the 
softening of an optical mode linearly coupled to a transverse acoustic phonon, 
the critical exponents retain their mean-field values. Disorder may lead to a va- 
riety of precursor effects and modified critical behaviour. Defects that locally 
soften the crystal may induce the phenomenon of local order parameter conden- 
sation. When the correlation length of the pure system exceeds the average defect 

— 1/3 

separation n D , a disorder-induced phase transition to a state with non-zero 
average order parameter can occur at a temperature T c (n-£,) well above the tran- 
sition temperature T® of the pure crystal. Near T®, the order-parameter curve, 
susceptibility, and specific heat appear rounded. For T < T c (no) the spatial in- 
homogeneity induces a static central peak with finite q width in the scattering 
cross section, accompanied by a dynamical component that is confined to the 
very vicinity of the disorder-induced phase transition. 
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1. Introduction 

Generally, displacive structural phase transformations (for reviews, see Bruce & 
Cowley 1981 and Rao & Rao 1987) can be divided into two different groups, 
namely distortive and elastic phase transitions. At distortive phase transitions 
some of the ions or molecular groups in the crystal's elementary cell are displaced 
with respect to each other. The corresponding displacement field serves as an ap- 
propriate order parameter for the transition, which in the case of a second-order 
transition is thus characterized by the softening of an optical phonon mode that 
becomes overdamped in the vicinity of the transition. Ionic displacements may 
lead to the appearance of a macroscopic polarization, i.e., to ferro- or antiferro- 
electric behaviour. Famous examples of this group of structural phase transitions 
are the perovskites with the high-temperature simple-cubic ABO3 crystal struc- 
ture; e.g.: the ferroelectric BaTi03, the antiferroelectric NaNb03, and SrTiO"3, 
which displays an antiferrodistortive transition accompanied by the softening of 
a zone-boundary phonon (see Liithi & Rehwald 1981). At these distortive trans- 
formations the elastic degrees of freedom constitute merely secondary variables, 
their interaction with the critical order parameter fluctuations being quadratic 
in the order parameter and linear in the strains. 

At an elastic phase transition (also called ferroelastic) , on the other hand, 
the crystallographic unit cell undergoes an elastic deformation, as is depicted in 
Fig. [I] for NaOH (for a review, see Wadhawan 1982). Thus the order parame- 
ter in this case is a certain linear combination of appropriate components of the 
strain tensor. In the majority of cases the crystal undergoes a shear deforma- 
tion, and the accompanying soft mode is the corresponding transverse acoustic 
phonon. KH2PO4 (KDP), KCN, and Fe30"4, for example, display a first-order 
elastic phase transformation, while in K-Na-tartrate, NbsSn, TeC-2, KH3(Se03)2, 
LaPsOi4, and DyV04 the transition is continuous (second-order), which im- 
plies that (at least) one of the sound velocities vanishes at the transition (Liithi 
& Rehwald 1981; Cummins 1982). As a consequence of the crystal anisotropy, 
the sound velocities depend on the direction of propagation, and thus the soft- 
phonon velocity does not vanish throughout the entire Brillouin zone, but only in 
one or several so-called soft sectors which may be one- or two-dimensional. E.g., 
for the martensitic transformation in the A-15 compounds NbsSn and the 
transverse acoustic phonons propagating along the face diagonals of the Brillouin 
zone with sound velocity c s = \J (en — cy^jlp soften, see Fig. I (Rehwald et al. 
1972); in KCN the elastic constant C44 vanishes at the transition, and the soft 
subspaces are two-dimensional (Haussiihl 1973; Knorr et al. 1985). We remark 
that the actual situation may be quite complex; e.g., in the ferroelectric KDP 
the z component of the electric dipole moment couples linearly to the £12 shear 
deformation (Brody & Cummins 1974). Consequently, although caused by the or- 
dering of hydrogen-bonding protons, the transition becomes ultimately an elastic 
phase transformation with vanishing shear modulus cqq. Similar cases, where the 
instability itself occurs for an optical mode that then drives the elastic transition, 
are LaPsOu (Fox et al. 1976) and h-BaTi03 (Yamaguchi et al. 1995), amongst 
others. Yet there are also cases of "pure" elastic phase transitions, where local 
fluctuations render the crystal unstable, e.g., in NaOH (Bleif et al. 1971). The 
other possibility, realized in the so-called isostructural elastic phase transforma- 
tions, is related to those elastic stability limits where all transverse and necessarily 
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all longitudinal sound velocities remain finite. For instance, if the bulk modulus 
of an isotropic or cubic elastic medium vanishes, only the macroscopic uniform 
dilatation and gradient modes soften, but none of the phonons (Khmel'nitskii 
1975; this kind of behaviour was observed in Ce, see Poniatovskii 1958). 

Figure |l] and Figure [2] near here 

In the present brief overview, we shall focus on the theory of continuous elas- 
tic phase transitions accompanied by the softening of an acoustic phonon. We 
shall rely on the Ginzburg-Landau free-energy functional that was established 
by Cowley (1976) and Folk et al. (1976 a, b), who also derived the Langevin 
equations of motion for the soft acoustic phonons (Folk et al. 1979). We are go- 
ing to investigate the statics and soft -phonon dynamics of second-order elastic 
phase transitions in pure systems in § 2 and § 3, respectively; for earlier reviews, 
see (Schwabl 1980) and (Schwabl 1985). As a consequence of the anisotropy of 
the elastic system, fluctuations are reduced and the (upper) critical dimension 
d c (m) = 2 + m/2 (for the case of a m-dimensional soft sector). This is smaller 
than d c = A for a (j) 4 model, describing either an Ising or Heisenberg magnet, 
or distortive structural transitions with short-range interactions. If the theory is 
phrased in terms of deformations, the Hamiltonian is seen to be equivalent to a 
spin system with long-range uniaxial dipolar interactions, which again explains 
why the corresponding critical dimension is lower (see §2). The case where the 
elastic transition is actually driven by the linear coupling of transverse acoustic 
modes to a softening optical phonon will be briefly discussed in § 4. 

The at the time surprising observation of extremely narrow central peaks in 
the scattering cross-section, near both distortive (Riste et al. 1971; Shapiro et al. 
1972) and elastic transitions (Shirane & Axe 1971) has prompted various theo- 
retical studies investigating the influence of lattice defects on the statics and dy- 
namics of structural phase transformations (e.g., Folk & Schwabl 1974; Halperin 
& Varma 1976; Hock & Thomas 1977; Schmidt & Schwabl 1977, 1978; Hock et 
al. 1979; Sasvari & Schwabl 1982; Weyrich & Siems 1984; Wiesen et al. 1987, 
1988). For reviews of the experimental facts, see (Miiller 1979) and (Fleury & 
Lyons 1982); some key theoretical results are collected in (Bruce & Cowley 1981). 
Recently, these studies, which are typically concerned with distortive transforma- 
tions and single-defect properties or conclusions derived by linear superposition 
thereof, were extended to elastic systems (Schwabl & Tauber 1991 a), and the 
methods developed for distortive crystals with a truly finite concentration of im- 
purities (Schwabl & Tauber 1991 b) were applied to the highly anisotropic elastic 
structural phase transitions (Bulenda et al. 1995). These developments and their 
consequences for a possible explanation of the central-peak phenomenon are the 
subject of § 5. In the concluding § 6 we summarize by contrasting elastic phase 
transitions with distortive structural transformations, and also discussing some 
aspects of the related first-order martensitic transformations. 

2. Static properties in pure systems 

(a) Ginzburg-Landau free energy functional 

In order to study elastic phase transitions, we expand the elastic free energy 
in terms of the strain tensor (Einstein's sum convention is employed here and 
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in the following), 



1 (2) 1 (2) / d \ ( d \ 

2 C iklm £ ik £ lm + ^""iklmrs \Q^T £ik ) \dx~ £lm ) ^ ' 

1 (3) 1 (4) 

+ ^ C iklmrs e ik^lm^rs + ~^ c iklmrsuv ^ik^lm^rs^uv + • • • 



d 3 x 



where the strain tensor is defined via the derivatives of the displacement field Uj , 

lk 2 \dxk dxi dxidxk) ' 

and the displacement field itj(x) itself can be expanded in terms of the Fourier- 
space normal coordinates Qk,A (A = 1,2,3), 

Ui (x) = ^=^e jkx e,(k,A)Q k , A . (2.3) 
V^ M k,A 

Here, e(k, A) denotes the polarization vector, M the mass and iV the number of 
the unit cells. 

The coefficients c (m) in Eq. 

are the isothermal elastic constants of order 



m. In addition to the terms contained in the standard expansion of the free 
energy in elasticity theory, we explicitly take into account contributions involving 
gradients of the strain tensor. Contributions of this form were introduced in so- 
called generalized elasticity theory, and can be derived from a lattice-dynamical 
model by taking the long-wavelength limit and retaining terms up to the fourth 
order in the wave vector k (see Krumhansl 1968; Mindlin 1968; Kunin 1968). 
These gradient terms are of course essential for the mathematical description of 
spatial fluctuations and short-range inhomogeneities. 



Figure 3 and Figure 4 near here 



The point -group symmetry now determines the number of independent elastic 

(2) 

constants c iklrn and, via the stability limits, the possible elastic phase transforma- 
tions (see Aubry & Pick 1971, Liakos & Saunders 1982). E.g., in the orthorhombic 
system one finds elastic phase transitions with vanishing C44, C55, and cqq (we ap- 
ply the Voigt notation). The corresponding wave vectors and polarization vectors 
are shown in Fig. |||a. In case of higher symmetries softening may occur in one- 
or two-dimensional sectors (if C44 — » 0) of Fourier space. The possible acoustic 
soft modes in tetragonal, cubic, and hexagonal crystals are illustrated in Figs. || 
and ||, respectively. Since hexagonal crystals are elastically isotropic with respect 
to the six-fold c axis, the transverse phonons with wave vector and polarization 
perpendicular to the c axis are degenerate and soften with cqq = (en — ci2)/2 
going to zero. 

If there are third-order invariants which lead to a cubic term in the soft phonon 
mode, the transition will be of first order. The third-order invariants allowed by 
symmetry were tabulated by Brugger (1965). E.g., in cubic systems with C44 — > 
there is a third-order term of the form £i2£23£i3; third-order terms are also 
present for c\\ — cyi — > in cubic and hexagonal systems. The theory described 
in the following is applicable, if third-order terms are either not present at all, 
or if they are sufficiently small in order that the ensuing first-order character of 
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the transition will only be noticeable in the immediate vicinity of the transition 
temperature T c , and the phase transformation can effectively be considered as 
continuous. The high-temperature phases and vanishing combination of elastic 
constants for the possible elastic phase transitions are listed in table |l|, along 
with the corresponding strain components [e 2 = (en — £22) / 'V%, 63 = (en + £22 — 
2e33)/v / 6] ) dimensionality m of the soft sectors, and the third-order invariants. 
For a listing of physical examples for these transformations, see (Folk et al. 1979; 
Liithi & Rehwald 1981; Cummins 1982). 



Table 1 near here 



A few remarks are in place here concerning the soft-mode spectrum and the 
char acteristic Hamiltonian or free-energy density following from the expansion 
(|2.lD, as to be investigated in the following § 2 6 and § 3. 

(i) The sound velocity vanishes only in one- or two-dimensional subspaces. In 
the vicinity of these directions the sound velocity is non-zero but small, and it 
is important for the theory to include all phonons with wavevectors k in a finite 
sector around the m-dimensional soft subspace. 

(ii) The gradient terms stemming from generalized elasticity theory prevent 
the sound frequency from vanishing throughout the entire Brillouin zone, and 
replace the linear by a quadratic dispersion precisely at the critical temperature 
T c . 

(iii) For the study of critical phenomena, we shall disregard non-critical modes 
and only retain the normal coordinate of the soft acoustic phonon. 

(iv) We shall discard odd anharmonic terms, thus restricting the applicability 
of the theory to situations where these are prohibited by symmetry; this applies 
to orthorhombic and tetragonal crystals, or to cases where the cubic terms are 
sufficiently small such that the phase transition can be regarded as of nearly 
second order. 



(b) Critical statics 

We are now ready to address the static critical behaviour near elastic phase 
transitions. In the framework of renormalization group theory, one considers a 
d-dimensional crystal with an m-dimensional soft subspace (in real systems, 
of course, d = 3 and m = 1 or m = 2). Accordingly, we decompose the d- 
dimensional wave vector k into its m-dimensional "soft" components q and the 
(d — m)-dimensional "stiff" components p; k = (q, p). The Hamiltonian following 
from Eqs. ( |2.1| ), fl2.3| ) and the discussion at the end of the preceding Section then 
reads in Fourier space (Folk et al. 1976 a, b) 

H[{Q*}] =lJd d k(rp 2 + p* + q 2 )\Q k \ 2 

+u J d d k l ...d d hv(k 1 ,... ,k 4 )Q kl Q k2 Q k3 Q k4 , (2.4) 

where r oc T — T c is assumed (in the following, we shall omit the polarization 
indices A). We will investigate the following two models: 

I : w(ki, k 2 , k 3 , k 4 ) = p\P2PWi <5(ki + k 2 + k 3 + k 4 ) , (2.5) 
II : u(ki,k 2 ,k3,k4) = (piP 2 )(P3P4)£(ki +k 2 + k 3 + k 4 ) . (2.6) 
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The characteristic features of the elastic Hamiltonian ( |2.4D are 

(i) the anisotropy in the harmonic part, as a consequence of which fluctuations 
in the "stiff" directions are suppressed, and 



(ii) the wave-vector dependence of the interaction (2.5) or Q2.6| ). 

We note that the statics of model I can be mapped onto the Hamiltonian of 
the uniaxial dipolar magnet using the transformation = pQk (Cowley 1976; 
Folk et al. 1977). 

Terms which are irrelevant for the critical behaviour (in the renormalization 
group sense) have already been omitted in the effective Hamiltonian ( |2.4| ) and 
interactions ( |2.5| ),( |2~6| ). As noted above, there is in general more than one soft 
sector; however, the respective interaction vanishes upon repeated application of 
the renormalization-group transformation. Hence we may consider the different 
soft sectors as independent and of the form ( |2.4| ). 

The renormalization-group transformation appropriate for the anisotropic 
Hamiltonian ( |2.4| ) consists, as usual, of two steps: 

(1) Eliminate wave vectors in the momentum shells b^ 1 < p < 1 for the soft, 
and b~ 2+v ^ 2 < q < 1 for the stiff sector, respectively. 

(2) Rescale according to p' = bp, q' = 6 2_ ' 7//2 q, and Q' k , = C _1 Qk- 

The elastic anisotropy is reflected in the different scaling for the soft and stiff 
sectors. This procedure yields a new effective Hamiltonian with coupling con- 
stants r' and u'. The decisive transformation is the one for the nonlinear coupling 

u, 

v! = b i+m ~ 2d \u - 36u 2 C(r)l , (2.7) 



where C(r) stems from the one-loop bubble diagram; from Eq. (2/7) we can imme- 
diately read off the (upper) critical dimension as a function of the dimensionality 
m of the soft sector (Folk et al. 1976 a, b) 

vn 

d c {m) = 2 + — . (2.8) 

For d > d c (m), the nonlinear coupling is irrelevant, and the system approaches 
a Gaussian fixed point u* = 0; its critical behaviour is thus governed by the 
classical critical exponents. Only for d < d c (m) is a non-trivial fixed point of 
the renormalization-group transformation approached, and the exponents assume 
non-classical values. 

The effective suppression of critical fluctuations now becomes evident. For one- 
dimensional soft sector the critical dimension is d c (l) = 5/2, and consequently 
elastic phase transitions in three dimensions are characterized by the classical 
critical exponents 

m = l: v = l,n = 0, 7 = 1, f3=-,a = 0,6 = 3 . (2.9) 

In the case of a two-dimensional soft sector one has d c {2) = 3, and one thus 
finds classical behaviour with logarithmic corrections in three dimensions. For 
instance, the static susceptibility (inverse elastic coefficient) and the specific heat 
are given by 

m = 2: x oc T^lm-Tp , C oc|lnT| rc , (2.10) 
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where r = \T — T c \/T c . The exponents r x and rc for models I and II are listed in 
table |2] (Folk et al. 1976 a, b). For a second-order elastic transition in d = 3 with 

m #S 2, the local fluctuations diverge: (it 2 ) oc | logr| for m = 2 and (u 2 ) oc r -1 / 2 
for m = 3; hence the Debye-Waller vanishes at T c . It is conceivable that the 
actual phase transformation in such systems will be of first order, not necessarily 
leading to the phase indicated by the soft mode. 

Elastic phase transitions with a two-dimensional soft sector occur, e.g., in KCN 
(Haussiihl 1973), NaCN (Rowe et al. 1975), and s-triazene (Smith & Rae 1978; 
Rae 1978); in most of the cases, however, the soft sector is one-dimensional. 



Table 2 near here 



There seems to be no isotropic elastic system that would show the onset of a 
shear instability, possibly as a precursor to melting. However, it is conceivable 
that the case m = d could apply to phase transitions in disordered systems 
like polymer gels, which on a sufficiently large length scale can be viewed as 
isotropic elastic media. For isotropic elastic transformations the upper critical 
dimension would be d c (3) = 7/2, and one would have non-classical behaviour in 
three dimensions. For elastic pha se transitions of layers deposited on (preferrably 
amorphous) substrates, Eq. ([2.8D predicts non-classical behaviour. 



3. Dynamic properties 

(a) Soft acoustic phonons 

In order to derive the equations of motion for the so ft a coustic phonons, we 
construct an effective Lagrangean from the free energy (|2.1| ) and a kinetic term, 
reexpressed in terms of the displacement fields via 

L[{u t }\ = J I {^^) 2 d^xdt - T[{ Ul }] ; (3.1) 

here, p denotes the mass density of the unit cell. Following nonlinear (generalized) 
elasticity theory the deterministic part of the equations of motion can now be 
obtained by applying the variational principle 5L[{ui}]/5ui(x,t) = 0. Degrees of 
freedom other than the displacement field will lead to damping and noise; thus 
our complete Langevin-type equation of motion for the phonon normal modes 
becomes (Folk et al. 1979) 

MQ k = - 6H ^ Qk}] -Mr k Q k + r k . (3.2) 

For acoustic phonons, the damping coefficient reads 

T k = Dp 2 +Dq 2 , (3.3) 

where the damping constants D and D for the soft and stiff sectors, respectively, 
are different in general. The stochastic force r k results from the non-critical de- 
grees of freedom; it has zero mean and its fluctuations are related to the damping 
coefficient by an Einstein relation 

(r k (i)Tv(i')> = 2T k k B T5(k + k')5(t-t') , (3.4) 



Phil. Trans. R. Soc. Lond. A (1996) 



8 



F. Schwabl and U.C. Tauber 



which guarantees that the equilibrium distribution is 



P[{A»Qk}] «exp 



(3.5) 



where denotes the canonical momentum conjugate to Qk- 

For T > T c , the dynamic phonon susceptibility in mean-field approximation 
follows immediately from the Fourier-transformed equation of motion ([T^) with 
( |3.3| ) , and the harmonic part of the Hamiltonian (^4|) , 

X(k, w) = [-Mu 2 - iMoj(Dp 2 + Z)g 2 ) + rp 2 + p 4 + <? 2 1 ~ . (3.6) 



The poles of %(k, u;) determine the phonon dispersion relation. Eq. (|3.6[) thus 
shows that as T — > T c (r = £ -2 — > 0), the dispersion of the soft acoustic modes 
(q = 0) turns from linear to quadratic; assuming that D stays constant, the soft 
mode finally becomes overdamped in the vicinity of T c . The dynamical phonon 
correlation function -D(k, to), which is directly measured in scattering experiments 
as part of the dynamic structure factor 5(k, ui) (see §5 6), is related to the dy- 
namic susceptibility via the (classical) fluctuation-dissipation theorem 

D(k,u,) = ^Im X (k,u;) ; (3.7) 

to 

it is depicted in Fig. || for q = and M 1 / 2 !) = 1.0, both in the hydrodynamic 
regime (p£i = 0.1) and in the critical region (p^2 = oo). 



Figure |5| near here 



Before turning to a more thorough analysis of the critical dynamics, we re- 
mark that, to be precise, the adiabatic elastic constants should enter the dynam- 
ical equations of motion, and not the isothermal elastic constants used in § 2. 
The general thermodynamic relation between adiabatic and isothermal elastic 
constants (of second order) is (Leibfried & Ludwig 1961) 

C iklm ~ C iklm = ^rPikfilm , (3-8) 

where (3^. = (daik/dT) e is the temperature-derivative of the stress tensor and 
C £ is the specific heat at constant strain. In certain symmetry directions the 
adia batic and isothermal sound velocities coincide; for instance, according to 
(3.8) cff — cff = c 1 ^ — cf 2 (switching back to Voigt's notation). However, in a 
cubic crystal off = cfi + T/3 2 /C £ , and this difference might be noticeable and 
important when data are actually fitted over a wide temperature range. 

(b) Critical dynamics 

In situations where the soft mode is a propagating hydrodynamic mode, con- 
siderable information can be gained already by simply using the hydrodynamic 
result in conjunction with dynamical scaling. Thus one finds for the sound velocity 

Cs = ^c~fp oc r 7/2 oc £-( 2 -"V 2 , (3.9) 

where p is the mass density, c tx r 7 denotes an appropriate combination of elastic 
constants, i.e., the inverse of the static critical susceptibility, £ oc t~ u is the 
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correlation length, and the scaling relation 7 = u(2 — rj) has been used. Therefore 
we expect for the dynamical critical exponent z = 2 — rj/2. 

Indeed, supplementing the static rescaling operations of § 2 b by uj' = b z tu, and 
carrying out the dynamic renormalization group program, to one-loop order the 
following recursion relations for the dynamical parameters M and D are found 
(Folk et al. 1979): 

M' = b i - r '- 2z M , (3.10) 
M'D' = b 2 ^^ z \MD + u 2 E(D, M) In bj . (3.11) 

Here, the explicit form for E(D, M) may be inferred from the one-loop integral, 
but is not important for what follows. At the Gaussian fixed point (u* = 0, hence 
77 = 0) both relations lead to 

z = 2 , (3.12) 

as expected. In the case of a one-dimensional soft sector, furthermore, the validity 
of the mean-field results implies that 

c s oc \T-T C \ 1 ' 2 , (3.13) 
while the damping coefficient is temperature-independent, 

L>oc|T-r c |° ; (3.14) 
consequently the sound attenuation coefficient becomes 

Duj 2 



O.s 



2c? 



oc \T -T c \~ 3/2 uj 2 . (3.15) 



For m = 1, the dynamical susceptibility is given by the mean-field expression 
(3.6), see Fig. |5[ In the case of a two-dimensional s oft se ctor, r = £ -2 has to be 



replaced by the inverse static susceptibility in Eq. ( |2.10| ). The power law ( |3~T3l) 
has been confirmed by many experiments; the prediction for the phonon damp- 
ing was verified both by Brillouin scattering (Errandonea 1981) and ultrasonic 
attenuation experiments (Garland et al. 1984). 

Although this is of no relevance in real crystals, it is interesting to note that 
dynamical scaling breaks down in the hypothetical isotropic elastic system for 
d < d c . The origin of this peculiar behaviour is the appearance of a dangerous 
irrelevant variable, namely X = M l l 2 D, whose fixed-point value diverges. Under 
these circumstances there is no definite dynamical exponent; in the immediate 
vicinity of T c one finds z = 2 + cry (as for the relaxational model A), while in the 
hydrodynamic region the exponent for the sound velocity is z = 2 — rj/2, but yet 
another value emerges for the damping coefficient (Folk et al. 1979). 



4. Linear coupling of an optical to an acoustic phonon 

As mentioned in the Introduction, there are frequently situations where the 
driving mechanism for the elastic instability is the softening of an optical phonon 
which couples linearly to a transverse acoustic mode, e.g. in KH2PO4 (KDP) and 
LaP50i4. Then one has to study the interaction of the acoustic and the opti- 
cal modes. The ensuing coupled system can be described alternatively by modes 
which are linear combinations of these phonons; the actual soft mode is primarily 
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acoustic and is governed by the characteristic anisotropic elastic Hamiltonian. 
Hence the resulting critical behaviour in three dimensions is classical, even if 
the optical modes have their origin in short-range interactions that would by 
themselves allow for non-classical critical fluctuations. In KDP, there is a lin- 
ear coupling of the polarization to the shear £12 (Brody & Cummins 1974). In 
LaPsOu, a Raman-active optical phonon coupled to £13 drives the transition 
(Fox et al. 1976). 

A simplified Hamiltonian containing all the essential features is (Schwabl 1980) 

#[{Qk},{Qk}] = ^ U(k) 2 |g k | 2 +w a (k) 2 |g k | 2 + 5 ^Q k g_ k 
1 k 

+ n E---E^ 1 ---9k 4 A(k 1 + ... + k 4 ) , (4.1) 

ki k 4 

where A(k) = 1 if k equals a reciprocal lattice vector and zero otherwise. Q k and 
g k are the normal coordinates of the acoustic and optical phonons, respectively, 
and their uncoupled dispersion relations are 

u (k) 2 = a'(T-T^/T^ + c 1 k 2 x + c 2 k 2 y + c 3 k 2 , (4.2) 

and 

w a (k) 2 = a\k 2 + a 2 k 2 + a 3 k 2 . (4.3) 
The third and fourth term in the Hamiltonian (|4.lD represent the linear coupling 



of the acoustic and the optical phonon, and the anharmonic interactions of the 
optical phonons, respectively. Note that: 

(i) The transition is driven by the optical phonon, Eq. ( |4.2| ). Due to the short- 
range interactions, the optical phonon by itself would display non-classical critical 
behaviour. 

(ii) In KDP there are (uniaxial) dipolar forces present, and uo(k) 2 thus contains 
terms oc k 2 /k 2 ; hence, without the coupling to the acoustic phonon, one would 
find logarithmic corrections to the classical exponents in three dimensions. 



(hi) Eqs. (4.1)— (4.3) refer to the vicinity of the k x axis. If th e aco ustic phonon 
couples to £12, there is also a contribution oc ikyQ^q^ in Eq. (4.1). 



The dynamics of the coupled system are described by the stochastic equations 
of motion 

ggKgkMQk}] r ( . v ^ (AA s 

Qk = 7 r, k g k + r 9k , (4.4) 

Qk = 77; r Q (k)Q k + r Qk . (4.5) 

The stochastic forces r q k and tq k are related via Einstein relations of the form 
( |3.4|) to the damping coefficients, which have the form lim k ^or g (k) 7^ and 
r<2(k) = D\k 2 + D%ky + D 3 k 2 . For small wave vectors, the modes which diag- 



onalize the harmonic part of the Langevin equations of motion (4.4), ( |4.5| ) are 
given by 

q k - a(k)Q k _ -a(k)g k + Q k 
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with a(k) = igk x /2uo(0) , and the corresponding eigenfrequencies are 

X ^"^ )2± M? ■ (4 ' 7) 

For small wave vectors k, the mode /2 k is essentially acoustic; its sound velocity 
in the k x direction, as obtained from A2, vanishes at the temperature 

T o 2 

r c = T c ° + -^ 7 . (4.8) 
In terms of the new variables /i/2k the Hamiltonian reads 

ff[{/i/2k}] = \ E [ A i l/ikl 2 + A 2 |/ 2k | 2 ] (4.9) 
z k 



\ - Aki + q(ki)/2ki /ik4 + «(k 4 )/ 2 k 4 , ( , , , , 



k 4 



Vl + I«(ki)| 2 v / l + l«(k 4 )| 2 



■1 



It is important to realiz e th at the terms involving /2 k ar e governed by the 
anisotropic Hamiltonian ( |2.4D with m = 1 and the soft direction k x . Since Ai 
remains finite at T c , /ik scales as /ik °c b 2 f2k an d is thus irrelevant for the 
critical behaviour. Therefore, the statics and dynamics of the phase transition at 
T c are described by mean-field critical exponents (Schwabl 1980), as previously 
suggested by a self-consistency argument (Villain 1970). 

Closing our discussion of pure crystals, we remark that interactions of the 
soft acoustic phonon with other acoustic modes have been disregarded here. For 
instance, terms of the form en x (shear) 2 are possible, which cause a coupling 
of the soft transverse mode to longitudinal phonons, similar in structure to the 
magnetostrictive interaction in compressible magnets (Wegner 1974; Bergman 
& Halperin 1976). Eliminating the non-critical longitudinal fluctuations leads 
to additional (negative) fourth-order couplings betweeen the soft modes. One 
may thus anticipate that in analogy to the magnetostrictive case, the transition 
becomes of first order for free boundaries, and the specific-heat exponent will be 
positive, a > 0. For a clamped crystal, on the other hand, the transition would 
remain continuous; however, clamping would interfere with the shear deformation 
which is characteristic of the elastic systems under consideration here. 



5. Disorder effects 

(a) Local order parameter condensation 

Real crystals contain defects of all sorts: point disorder, dislocations, grain 
boundaries, etc. The more complicated the unit cell, the more likely crystal 
growth will lead to imperfections. Some of these defects will have a pronounced 
impact on the phase transition. We thus turn to the investigation of the influ- 
ence of disorder on the statics and dynamics of elastic phase transitions. More 
specifically, we shall be interested in a certain type of inhomogeneities, which 
have the effect of locally changing the elastic constants, and thus softening the 
crystal. In this Section, we shall illustrate some local features by discussing an 
effectively one-dimensional model with a single impurity (Schmidt & Schwabl 
1977, 1978; Schwabl & Tauber 1991 a). As a matter of fact, three-dimensional 
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systems behave qualitatively similar, when we restrict ourselves to the soft sectors 
in momentum space. Understanding these single-defect properties facilitates the 
interpretation of the results for systems with a truly finite disorder concentration, 
to be addressed in §56 (Schwabl & Tauber 1991 b; Bulenda et al. 1995). 

The one-dimensional elastic model free energy reads (Schmidt & Schwabl 1978) 

where the harmonic elastic constant vanishes at the transition temperature of 
the pure crystal, a = a'(T — T c °), and the Ginzburg-Landau parameters b and 
c are taken to be positive constants. The disorder influence is described by the 
short-range defect potential U(x), which can be approximated by a delta func- 
tion, U(x) = Uq6(x), if both the correlation length and the wavelengths of the 
characteristic fluctuations are small compared to the typical defect size, which 
are assumed to be comparable to the lattice constant clq. An attractive potential 
with positive Uq means that the transition temperature is locally increased. The 
probability distribution V[e\ for a configuration e{x) is given by 

V[e]otexp(-F[£]/k B T) ; (5.2) 

in the spirit of Ginzburg-Landau theory one determines the most probable state 
e(x) from the stationarity condition 5J-\e\/ 5e{x) = 0; this one-dimensional con- 
figuration then serves as an approximation for the real three-dimensional order 
parameter (for a review of the application of Landau theory to structural phase 
transitions, see, e.g., Salje, 1992). 



Inserting (5.2) yields the nonlinear second-order differential equation 

cs(x)" = [a — U(x)]e(x) + be(x) 3 (5.3) 

for the stationary state, which of course always has the constant solution £q = 
(note, however, that the constant solutions for T < T c °, e± = ± y/\a\/b describing 
the two possible orientations of the ordered state, are no longer allowed if Uq ^ 0). 
Due to the defect influence this homogeneous state may become unstable at a 
certain temperature T l c > T c °, and instead a localized order parameter conden- 
sate (cluster) near the defect may form (Schmidt & Schwabl 1978); T l c is called 
local transition temperature, although it does not define a proper phase tran- 
sition. Linear stability analysis shows that in the one-dimensional model local 
condensation occurs at 

T l c = T° c +^- , (5.4) 
Aa'c 

for any positive defect strength Uq. For T < T l c one finds the following stable 
cluster configurations (see Fig. ||) in the interval T® ^ T ^ Tj, 



£>{x) -sinH\x\/C >+P> ) ' (5 - 5) 
where £> = \Jcja and p> = arcoth(f7o/2y / ac); and in the temperature range 



T <T C 



o 



£<(*) = ± v /|a|/6coth(|x|/£<+ /0< ) , (5.6) 
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with £< = y/2c/\a\ and p< = (1/2) arsinh(^/8|a|c/J7o)- Note that the widths of 
both configurations are given by the correlation lengths of the pure system £> 
and £<, respectively. We remark that in d > 1 dimensions, in general a certain 
minimum disorder potential strength must be exceeded for local condensation 
to occur (Schmidt & Schwabl 1977; Schwabl & Tauber 1991 6); but the ensuing 
order parameter cluster still decays exponentially oc e~ r ^ for large r S> £. 



Figure |6| near here 



One can now look for traces of the local order parameter condensation in the 
spectrum of the soft phonons. The one-dimensional model ( |5.1[ ) applies for both 
distortive and elastic transitions, and hence the static single-defect properties 
are completely analogous. However, their respective dynamics are very different: 
In the case of distortive transitions, one may obtain the spectrum of the optical 
phonons for T ^ T l c by linearizing Eq. ( |5.3| ) about the stationary solution Eq = 0; 
as a result of the "attractive" defect, one additional localized mode emerges below 
the continuum of propagating scattering states, and its frequency vanishes pre- 
cisely at the local transition temperature Q5.4|). This softening of a localized defect 
excitation can thus be viewed as a dynamical precursor for the local condensation 
in distortive systems (Schmidt & Schwabl 1978). 

Similar localized acoustic phonon modes have been hypothesized for elastic sys- 
tems, namely for the related first-order martensitic transformations (Cla pp 1973, 
1979, 1981). When supplemented w ith a dynamical equation of the form (|3.2| ), the 
simplified one-dimensional model (5.1) permits an analytic solution. Linearizing 
the equation of motion yields a fourth-order differential equation for the soft 
acoustic modes, the eigenstates of which are linear combinations of propogating 

(oc e ±ikx ) and exponentially localized (oc > x "j contributions. Employing 
the appropriate boundary conditions for the delta-function defect, one finds that 
there appears no additional mode in the elastic system; any dispersion-free lo- 
calized mode would be unstable towards decay into continuum states. Instead, 
each of the propagating modes contains a localized vibrational contribution in 
the defect vicinity, which for the long-wavelength modes (k — ► 0) condenses as 
T — ► Tj, in the sense that an incoming phonon is not transmitted through the 
defect region, but the local order parameter condensate is formed (Schwabl & 
Tauber 1991 a). Accordingly, the delay time of an incoming wave packet diverges 
oc (T — T l c )~ l as T l c is approached, as does the relaxation time of a stress-induced 
cluster above the local transition. 



(b) Disorder-induced phase transition and central peak 

By appropriately averaging the single-defect solutions, one may already infer 
some properties of a crystal with Nd S> 1 defects, the concentration n-Q = N^/N 
of which is to remain finite in the thermodynamic limit (Schwabl & Tauber 1991 
a). Instead, we proceed to study such a system in d space dimensions using a 
more direct approach (Schwabl & Tauber 1991 b; Bulenda et al. 1995). It is now 
more convenient to use a lattice representation, and we therefore define the short- 
range quenched impurity potential produced by Nj) point defects, located on the 
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randomly selected lattice sites in as 



7V D 
in=l 



with Fourier transform (note that the system is not translation-invariant) 



1 



N 



-i(kxj— k'xj) 



(5.7) 



(5i 



As in the previous Section, we require that A = Uq/ciq > such that the defects 
locally enhance the transition temperature. One can then write the harmonic 
part of the elastic Hamiltonian (|2,4j), supplemented by the disorder potential, as 

H[{Q*}] = [d d k[ d d k'\ \(ap 2 + aq 2 + c/)<W - (kk> kk >l Q k Q- k > + 0(C&) , 



(5.9) 

where a = a' (T - T c °), and a,b > are constants. In addition to the thermal 
average, for this system with random quenched disorder each physical quantity 
has to be averaged over all possible defect configurations; the formal definition 
of this quenched disorder average is 



((•••)) 



n 

i=i 



1 

N 



N 

E 

i Dj =l 



(5.10) 



The soft-phonon dynamics is given by the Langevin equation ( |3.2|) , with 
Eqs. ( |3.3| ) and (3.4). In the pure system, the ("free") phonon propagator reads 
[see Eq. (O)] 



-Muj 2 - iMuj{Dp 2 + Dq 2 ) + ap 2 + aq 2 + cp 4 



(5.11) 



Introducing an external field /i k in the equation of motion, inverting its Fourier 
transform, taking the thermal average, and differentiating it with respect to 
one arrives at the following recursion relation for the "full" dynamic response 
function in the high-temperature phase, 

X(k,k» = Xo (k,o;) ( 5 kk , + xo(k,a;)^(kk")^ kk »x(k /, ,k', W ) ; (5.12) 

k " 

system atica l iteration of Eq. ( |5.12| ), and then performing the configurational av- 
erage ( |5,10D finally yields the translationally invariant response function x(k, u) 
(for a diagrammatic representation, see Schwabl & Tauber 1991 b). Collecting all 
contributions that are linear in the defect concentration no (single-site approxi- 
mation), one eventually finds (Bulenda et al. 1995) 



X(k,w) 



Xo(k,w) 1 



Xnu(p 2 + q 2 



1 - A(a /2vr) rf J A (p 2 + q 2 ) X o(K co)d m pd d - 



(5.13) 

where the momentum cutoff A ~ 2-7r/ao has been introduced. 

The result ( |5,13| ) implies that as a consequence of the coupling to the softening 
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defects, the entire system may become unstable towards a new low-temperature 
phase with finite average order parameter at a certain temperature T c (nrj), which 
is to be determined from the condition 



lim 



k^ 2 x(k,u = OH 1 =0 . (5.14) 



As for distortive transitions (Schwabl & Tauber 1991 b), for d > 1 a certain 
minimum defect strength is required for this instability to occur, and T c (no) 
is bounded below by the local transition temperature T l c of § 5 a , which can be 
considerably higher than T®, the transition temperature of the pure system. 

In the response and correlation functions, the singularity at T c (ns) appears as 
a dynamical central peak. In view of the single-defect properties of the previous 
Section, one may understand this peak as stemming from overlapping localized 
vibrations at the point defects which "condense" at T c (no), thereby forming a 
spatially inhomogeneous order parameter configuration. In Fig. [7] the dynamical 
correlation function ( |3.7| ) is depicted for the case of an elastic system with a 
one-dimensional soft sector; the numerical values used for the Ginzburg-Landau 
parameters there are appropriate for NbsSn (Tj? = 45K), and the defect strength 
A was adjusted arbitrarily in order that T c (u-d) = 65K (Bulenda et al. 1995). As 
T — > T c {n-£>) from above, the soft phonon peak at finite energy is shifted to lower 
frequencies, and a very sharp and distinctive dynamical central peak emerges 
[Fig. |7|(a)]. Fig. |7](b) shows the dependence on the angle 9 between the external 
wave vector k and the soft sector; both graphs demonstrate that the dynamical 
central peak is confined to temperatures very close to the defect-induced phase 
transition at T c (nrj), and to wave vectors within the soft sector, reflecting the 
fact that wave vectors in the stiff directions do not probe the critical properties 
of the crystal. 



Figure [7| near here 



The properties of the ensuing inhomogeneous low-temperature phase may be 
studied using a self-consistent mean-field approach (Schwabl & Tauber 1991 b; 
Bulenda et al. 1995). The starting point is the following discrete version of the 
nonlinear Ginzburg-Landau free ener gy fo r a single-component order parameter 



£i in d dimensions [cf. Eqs. (2.1) and (5.1)1, 



^ N ^ N N D ^ N N 

Ffo] = 2 £ i G 0ij £ i ~ 2 X) £ i 6 hiD + 4 S e < - £ hi£i ' ( 5 ' 15 ) 
i,j=l i=l *d=1 i=l i=l 

where a stress term oc hi has been introduced, and the static propagator Gq^ is 
defined by its Fourier transform 

G ^ = ap2 + a £ + CP * . (5.16) 

In the framework of the Ginzburg-Landau approximation, i.e., neglecting order 
parameter fluctuations, the stationarity condition becomes (with constant exter- 
nal stress hi = h) 

G olj £ j ~ X J2 £i6i ' i o +bef = h . (5.17) 

3 *D 
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The solution of Eq. ( |5.17D , with its combined nonlinearity and randomness, 
in general poses a formidable problem. Therefore an additional approximation is 
employed, namely the following ansatz for the thermodynamical average of the 
order parameter is used, 

ei = A + Bj2$i,in > ( 5 - 18 ) 



i.e.: the order parameter at each lattice point i is assumed to be the sum of a 
homogeneous background A and an additional contribution B, if there is a defect 
at site i, thus enhancing the total value of the order parameter to A + B (see 
Fig. |8|). Thus we use the approximation that at all defect sites the order param- 
eter points in the same direction, and in addition neglect the spatial variation 
of the order parameter near the defects. This seemingly rather crude ansatz al- 
ready contains the possible relevant modifications which can be caused by the 
impurities, namely (i) an enhancement of the spatially averaged order parameter 
(corresponding to the parameter A), and (ii) the ensuing "screening" of the defect 
potential (described by B). 



Figure |8| near here 



Inserting the ansatz ( p. 18 ) into the stationarity equation ( 5.17| ) yields 



£k = hS K0 G (k) + G (k) ^2 <Akk'£k' , 

k' 

where a renormalized propagator 

Go(k)- 1 = G (k)- 1 + 6,4 2 
and screened defect potential (ft^i with weakened strength 

A = A - b \(A + B) 2 - A 2 



(5.19) 



(5.20) 



(5.21) 



have been introduced. As E q. fl5.12j ) for the dynamics in the high-temperature 
phase, the recursion relation ( |5.19j ) can be iterated and then the quenched disorder 
average performed, with the single-site approximation result 



m) 



bA 2 



An D 



(5.22) 



l-A(a /27r)% A G (k)d^J 

on the other hand, immediate averaging of Eq. (|5T7l ) gives 

(a + bA 2 )(A + n D B) -\n D (A + B) = h . (5.23) 

Eqs. ( 5.22 ) and ( |5.23 ) constitute two coupled nonlinear equations that uniquely 
determine the mean order parameter; assuringly, for h — ► they yield non-zero 
self-consistent solutions for {(s)} precisely when T < T c (ni)), the transition point 
determined from the instability of the high-temperature phase. Fig. |9] shows that 
the order parameter of the disordered crystal as function of T looks similar to the 
corresponding curve for the pure system, with the singularity at T® being smeared 
out by the defects (in the graph, the reduced temperature t = (T — T c )/T c ° is 
used). The order parameter sets in continuously at T c (nr>), with the usual mean- 
field exponent = 1/2, remains minute in the range T c (?t-d) > T > T®, and starts 
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to grow to larger values only near T®. Thus the transition temperature of the pure 
system remains an important parameter even in the perturbed system, although 
the true phase transition takes place at T c (nrj), which, however, may in fact 
be hardly noticeable in experiments. Again, the results for a three-dimensional 
system with one one-dimensional soft sector are depicted, but the qualitative 
features are essentially the same in the cases of a two- (or even three-) dimensional 
soft sector. 



Figure |9| near here 



Inserting the self-consistent solutions f or th e mean order parameter as func- 
tion of temperature into the free energy ( 5.15| ), the specific heat in Landau ap- 
proximation is readily determined by taking the derivatives with respect to T, 
C v = —T(d 2 J : [e]/dT 2 )v; see Fig. [l(| Obviously, the discontinuity at T<? has been 
smeared out, in place of which a tiny jump emerges at T c (jid). Note that there is 
a very distinct maximum of the specific heat near T^, while the extremely minute 
jump at T c (ni)) might not be experimentally detectable at all. 



Figure 10| near here 



The dynamical phonon response and correlation functions may similarly be 
calculated following the above procedu res; using the ansatz ( |5.18| ) leads to a 
recursion relation precisely of the form ( |5.12 ), where modified parameters 

a^a + 3bA 2 , a a + 3bA 2 , A -> A - 3bB(2A + B) (5.24) 

are to be inserted. One finds that the dynamical central peak in the phonon cor- 
relation function rapidly disappears as the temperature is lowered below T c (nj}). 

However, there is a very interesting additional static contribution to the dy- 
namic structure factor, caused by the defect-induced spatial inhomogeneity of 
the order parameter. The dynamic structure factor, as measured in scattering 
experiments, is defined as 



S(k,u) 



„iu>t 



\ l,J=l 



e -ik[a 4 +Ui(t)] e ik[aj+Uj(0)] W 



(5.25) 



where aj denote the Bravais lattice sites, and Uj the displacements from these 
equilibrium positions. In order to evaluate ((S(k,uj))), we perform a cumulant 
expansio n for the combined thermal and configurational averages of the exponen- 
tials in ( 5.251 ) and expand to second order, decomposing the displacement fields 
Uj(i) into a static part Uj and a fluctuating contribution iii(t). Eventually one 
arrives at (Schwabl & Tauber 1991 b; Bulenda et al. 1995) 



((S(k,u>))) =2tt5{uj) 



g ap 



-2W 



+ 



J^k^D^ik, 

a/3 



-2W 



(5.26) 



where the first term represents the elastic Bragg scattering peaks occuring at the 
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reciprocal lattice vectors g, W is the Debye-Waller factor, and 

S f(V = ^E^"""^ " «^»«^») (5-27) 

yields an additional static contribution to the structure factor arising from elastic 
scattering from random variations of the local order parameter (Huang scatter- 
ing); finally, D Q ^(k,w) is the dynamical phonon correlation function ( |3.7| ) and 
( |5.13| ), discussed before, describing inelastic scattering processes. 

In the framework of the above self-consistent mean-field theory and single-site 
approximation, the Huang scattering amplitude is readily obtained observing that 
ek = kiik, and using iterations of Eq. ( 5.19Q , with the result 



^ (k ) = ^A 2 (( £ -))WG o( k)2 = (a + hA , ? m? k ^G^f , 

1 - \(a /2ir) d J A G (k>)d d k'] n ° 

(5.28) 

where in the final step the condition ( 5.22j ) (for h — > 0) was used. Note that 



the renormalized static susceptibility, which determines this additional elastic 
scattering amplitude, does not diverge at T c °; furthermore, this static central peak 
is characterized by a finite width in momentum space, 7 = \/(a + bA 2 )/c (if k lies 
in the soft sector), which vanishes at the defect-induced tr ansit ion temperature 
T c (n\y). Fig. [ll] depicts the Huang scattering amplitude (|5.28[) as function of 
temperature for different wave vectors k = p in the soft sector. The static central 
peak emerges at T c (n-£>), and grows to considerable intensity near T®, but becomes 
less prominent and widens in k space for T < T®, as the order parameter slowly 
approaches the spatially homogeneous configuration of the undisturbed system. 
The dynamical central peak found in the phonon correlation function near T c (nj)) 
may be viewed as its dynamical precursor in the critical region, as well as the 
static one due to the new Bragg peaks which of course persist through the entire 
low-temperature phase. 



Figure 11 near here 



At this point we have to comment on the applicability of our mean-field ap- 
proach. Although we have explained in § 2 that critical fluctuations only play 
a minor role for elastic phase transitions, the situation here is somewhat differ- 
ent, for we now have to worry about local fluctuations in the orientations of the 



disorder-induced order parameter clusters (5.5), which were entirely neglected in 
the above treatment. The mean-field approximation suggests that as soon these 
local condensates form at T c {nj)) ~ T l c , a coherent order parameter with non— zero 
mean emerges, independent of space dimension d. More realistically, one would 
expect that while at T l c local condensation at the defects takes place, yet the 
orientations of the different clusters still fluctuate considerably, and they pos- 
sibly form a true new ground state only at a lower temperature T or( j. A crude 
estimate of this true ordering temperature is obtained by the criterion that the 
correlation length should at least be of the order of the mean defect separation 

oc ^rj 1 ^ for collective behaviour of the condensates to occur. A more refined ar- 
gument considers the different cluster orientations as effectively Ising-like degrees 
of freedom, and determines T or( j = J/k&, where J is the free-energy difference 
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between configurations with parallel and opposite cluster orientations (Bulenda 
et al. 1995). In d = 3 dimensions, the resulting "true" phase transition temper- 
ature is considerably lower than T c (71d) for identical disorder concentration, but 
may still be well above T®; thus, at least for T $J T or( j a spatially inhomogeneous 
order parameter appears, which produces a static central peak with finite width 
in momentum space for temperatures above the transition temperature of the 
pure system, where the average order parameter is still very small. In addition, 
the time scale of the cluster orientation fluctuations will diverge oc (T — T orc i) -1 
upon approaching T or( j, which in experiment would eventually render these in- 
distinguishable form static inhomogeneities. Qualitatively, at least, the situation 
will then appear as shown in Figs. f|-[ll], with T c {n^) replaced by T or d (or even a 
somewhat higher temperature, depending on the experimental frequency resolu- 
tion). 

(c) Disorder and critical properties 

Another important issue is the question if disorder may change the critical 
properties near elastic phase transitions, i.e., if the static and dynamical critical 
exponents are determined by a new renormalization-group fixed point. We briefly 
discuss the two relevant cases of (i) random-fields, and (ii) random-T c disorder, 
and their possible influence in the critical region very close to T c (Morgenstern 
1988). 

Quenched random fields couple linearly to the order parameter, and are taken 
to have vanishing average (/ik) and second moment 

(/ikM = A£T e 5(k - k') , (5.29) 

where = for short-range disorder correlation, while < < m for long-range 
correlated defects. Analyzing the renormalization-group recursion relations, one 
finds that these random fields may alter the pure result for the (upper) critical 
dimension ( |2.8j ) to 

m 

4(m) = 3 + - + - , (5.30) 

where m denotes the dimension of the soft sector (Morgenstern 1988). One there- 
fore expects non-classical critical exponents in three-dimensional crystals un- 
dergoing an elastic phase transformation under the influence of random strain 
fields. 

For disorder coupling quadratically to the order parameter, i.e., of the random- 
Tc type, on the other hand, a variant of the Harris criterion may be formulated. 
Again introducing possible long-range correlations 

(Sr^Sr^) = Wk~ e 6Qz - k') , (5.31) 

one finds that the second-order elastic phase transitions remains to be governed 
by the pure critical exponents, provided that 

a + z^0<O , (5.32) 

where a and v denote the specific heat and correlation length critical exponents, 
respectively (Morgenstern 1988). For the case of point defects (0 = 0) this means 
that new randomness-induced critical exponents emerge if a > 0, which is thus 
only possible for Ising-type systems. 
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6. Summary and comparison with related materials 

In conclusion, we summarize the above results and contrast them with related 
displacive structural phase transformations, namely (1) second-order distortive 
transitions, and (2) first-order martensitic transitions. 

For second-order distortive structural phase transitions, the order parameter 
is the displacement field corresponding to a soft optical phonon, while for elastic 
transformations the relevant collective variable is a certain combination of strain 
tensor components, and the soft mode is typically a transverse acoustic phonon. 
This means that while in the distortive case at T c one of the optical modes soft- 
ens at a specific k vector, an elastic instability implies the vanishing of the sound 
velocity along the soft sector in the Brillouin zone. As we have seen in § 2, as a 
consequence of the ensuing strong anisotropy critical fluctuations are suppressed, 
and continuous elastic phase transitions in three dimensions are governed by the 
classical critical exponents (for a one-dimensional soft sector) with at most log- 
arithmic corrections to the scaling functions (for a two-dimensional soft sector). 
For distortive structural transitions, on the other hand, critical fluctuations can 
be quite prominent and have important effects on the values of the critical ex- 
ponents (see, e.g., Bruce & Cowley 1981). Accordingly, the dynamical critical 
exponent is z = 2 for elastic transitions (with m < d), as in mean-field theory 
(§3), in contrast to the distortive case where typically z = 2 — crj with rj > 0. 
Elastic instabilities may actually be driven by the softening of an optical mode 
which is linearly coupled to the transverse acoustic phonons; but even in that 
situation the ensuing critical behaviour turns out to be classical; nevertheless, an 
anomaly is found in the sound attenuation coefficient (§4). 

We remark that the anisotropic Hamiltonian ( p.4| ) applies also to spin reori- 
entation transitions (Hornreich & Shtrikman 1976), to the phason instability at 
T = 49K in TTF-TCNQ (Bak 1976), and to the transition from smectic-A to 
smectic-C in a magnetic field (Hornreich & Shtrikman 1977). Anisotropic inter- 
actions are also present at Lifshitz points (Hornreich et al. 1975). 

The phenomenon of local order parameter condensation, induced by defects 
which locally increase the transition temperature, is qualitatively similar for both 
distortive and elastic phase transitions. However, dynamically the mechanism of 
condensation is quite different; while in the distortive case there appear localized 
modes below the continuum, which then soften at the local transition temperature 
Tg, there appears no such additional localized mode in the elastic case. Instead, 
the disorder leads to localized vibrational components in the propagating modes, 
which then "condense" at T l c to form the defect-induced order parameter clusters 
(§5 a). As was explained in §5 6, in a system with finite disorder concentration 
nj) a true defect-induced phase transition can occur at T c (nj)), resulting in a low- 
temperature state characterized by an initially strong spatial inhomgeneity, which 
finally becomes smoothened out at low temperatures, where the disorder influence 
is suppressed (§5 6). A dynamical precursor for this singularity at T c (no) (or 
T or d) is the appearance of a dynamical central peak in the critical region, which 
only in the case of distortive transitions can be traced back to the softening of 
a separate phonon impurity band. In the low-temperature phase, besides the 
new Bragg peaks an additional elastic Huang scattering component emerges as 
a consequence of the spatial variations of the local order parameter in either 
case. Near the transition temperature T® of the pure system, thermodynamic 
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quantities like the order parameter susceptibility and the specific heat appear 
characteristically rounded. Both these features have been observed in a variety 
of experiments. 

As was discussed in § 2, elastic instabilites often lead to first-order phase transi- 
tions, as is the case for the interesting class of martensitic transformations. These 
materials display a considerably more complex behaviour than the second-order 
transitions we have discussed here (the precursor effects in martensitic materials 
are reviewed in Finlayson, 1983). However, a simplifying elastic Ginzbur g-La ndau 
free energy has been proposed (Falk 1980, 1983), which is of the form (|2.1|) , but 
with a negative fourth-order term and including a sixth-order term in order to 
describe a first-order transition. Recently, an explanation for the experimentally 
observed "tweed" pattern in microscopic images of martensitic materials for tem- 
peratures far above the transition temperature has been suggested on a similar 
basis (Kartha et al. 1991, Sethna et al. 1992, Kartha et al. 1995). Namely, such 
pseudo-periodic lattice deformations were found to emerge in an elastic model 
including defects that locally modify the transition temperature in a random man- 
ner. This intrinsic compositional disorder was found to conspire with the natural 
geometric constraints of the lattice to form a frustrated, glassy "tweed" phase; 
and the random order parameter orientations then produce a static central peak 
with finite width in momentum space. We finally note that an interesting and 
wide field not covered in the present brief review concerns orientational glasses 
found in mixed-crystal solids (Hochli et al. 1990). 

Solid state physics in recent years has gone beyond the mere study of the 
materials offered by our surrounding nature, but more and more aims at the 
creation of new materials with taylored properties. As in semiconductor physics 
and magnetism, a similar development is about to take place in ferroelastic and 
ferroelectric substances. Besides of in-depth studies of disorder, the investigation 
of artificial multilayers and other composite structures will be of importance in 
the future. 
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High-temperature Vanishing combination Strain 
phase of elastic constants 



Third-order 
invariants 



Orthorhombic 



Tetragonal II 
Tetragonal I 



Cubic II 

Cubic I 
Hexagonal II 

Hexagonal I 



C44 
C55 

C66 
C44 
C44 
C66 
Cll — C12 
C44 

en - C12 

C44 
en — C12 
C44 

C66 = |(cil — C12) 
C44 

C66 = |(cil — C12) 



£23 
£13 
£12 
£23, £13 
£23, £13 
£12 
£11 — £22 
£12, £13, £23 

e3,e2 

£12, £13, £23 
e3,e2 
£23, £13 
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£23, £13 
£12, £11 — £22 



1 
1 
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1 + 2 
1 + 2 

1 

1 

2 

1 

2 
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1 + 2 
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1 + 2 
2 



£23 £13 £12 
e 3 (el-3e|), 
e 2 (e| - 3e^) 
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e 3 (ei-3ei) 



(£11 
£12 (£11 - 
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(£ll - £22) 
£l2(£ll 



£22) 



Table 2. Logarithmic corrections. 



Model 


r x 


rc 
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1/3 


1/3 


II 


4/9 


1/9 



Figure 1. Temperature dependence of the lattice parameters (a) and the monoclinic angle (5 

(b) in NaOH (from Bleif et al. 1971). 

Figure 2. Elastic moduli en, C12, and B of NbsSn vs. temperature (from Rehwald et al. 1972). 

Figure 3. (a) Wave vectors and polarization vectors of soft transverse phonons in orthorhombic 

crystals; the sound velocities c, are: {cn/ p) 1 ^ 2 (dashed), (css/p) 1 ^ 2 (solid), (cee/p) 1 ^ 2 (dotted), 
(b) and (c) Wave vectors and polarization vectors of soft elastic modes in tetragonal crystals. 

Figure 4. (a) and (b) Elastic soft modes in cubic crystals, with (a) wave vector in one of the 
base planes and polarization perpendicular to it, c s = (044/ p) 1 ^ 2 , and (b) wave vector in one of 
the face diagonals, c s = [(en — ci2)/2p] 1 ^ 2 ; (c) Soft transverse phonons in the hexagonal system, 
with sound velocities (C44/ p) 1 ^ 2 (dashed) and [(en — Ci2)/2p] x ^ 2 (solid). 
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Figure 5. Dynamical correlation function for the soft acoustic phonon vs. uiM 1 ^ 2 for q = 0, 
M 1//2 D = 1.0, and two values of the correlation length £; solid line: p£i = 0.1 (hydrodynamic 
regime), dashed line: = oo (critical region). 

Figure 6. Local condensation: Scetch of the order parameter profiles in a single-defect system 
in the temperature ranges T < T° , T° ^ T ^ T l c , and T > T l c , respectively. 

Figure 7. Dynamical phonon correlation function (with model parameters appropriate for 
NbsSn): (a) for different temperatures T and fixed angle 9 = 0, and (b) for fixed temperature 
T = 65. OIK with different angles 6 between external wave vector and soft sector (T c (no) = 65K). 

Figure 8. Schematic illustration of the order parameter configuration according to the 

approximations used here. 

Figure 9. Average order parameter ((e)) vs. reduced temperature t = (T — T c )/T c (a) 
(parameters appropriate for NbsSn); the region near t c (nn) = 0.444 is magnified in (b). 

Figure 10. Specific heat C v vs. reduced temperature t (a) (parameters appropriate for NbsSn); 
the region near t c (n-o) = 0.444, where the discontinuity occurs, is magnified in (b). 

Figure 11. Huang scattering intensity S c (k) vs. reduced temperature t for different wave 
numbers k = £2 3 / 2 7r/a ; the wave vector k lies in the soft sector. 
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